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Abstract. The curvature effect on a quantum dot with impurity is investi- 
gated. The model is considered on the Lobachevsky plane. The confinement 
and impurity potentials are chosen so that the model is explicitly solvable. 
The Green function as well as the Krein Q-function are computed. 

Keywords, quantum dot, Lobachevsky plane, point interaction, spectrum. 



1. Introduction 

Physically, quantum dots are nanostructures with a charge carriers confinement 
in all space directions. They have an atom-like energy spectrum which can be 
modified by adjusting geometric parameters of the dots as well as by the presence 
of an impurity. Thus the study of these dependencies may be of interest from the 
point of view of the nanoscopic physics. 

A detailed analysis of three-dimensional quantum dots with a short-range im- 
purity in the Euclidean space can be found in lj . Therein, the harmonic oscillator 
potential was used to introduce the confinement, and the impurity was modeled by 
a point interaction (5-potential) . The starting point of the analysis was derivation 
of a formula for the Green function of the unperturbed Hamiltonian (i.e., in the 
impurity free case), and application of the Krein resolvent formula jointly with the 
notion of the Krein Q-function. 

The current paper is devoted to a similar model in the hyperbolic plane. 
The nontrivial hyperbolic geometry attracts regularly attention, and its influence 
on the properties of quantum-mechanical systems has been studied on various 
models (see, for example, [2][3j[4]). Here we make use of the same method as in 
[T] to investigate a quantum dot with impurity in the Lobachevsky plane. We will 
introduce an appropriate Hamiltonian in a manner quite analogous to that of [T] 
and derive an explicit formula for the corresponding Green function. In this sense, 
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our model is solvable, and thus its properties may be of interest also from the 
mathematical point of view. 

During the computations to follow, the spheroidal functions appear naturally. 
Unfortunately, the notation in the literature concerned with this type of special 
functions is not yet uniform (see, e.g., [5] and [5]). This is why we supply, for the 
reader's convenience, a short appendix comprising basic definitions and results 
related to spheroidal functions which are necessary for our approach. 

2. A quantum dot with impurity in the Lobachevsky plane 
2.1. The model 

Denote by (g,<fi), < g < oo, < 4> < 2ir, the geodesic polar coordinates on the 
Lobachevsky plane. Then the metric tensor is diagonal and reads 

(gij) = diag^l, a 2 sinh 2 

where a, < a < oo, denotes the so called curvature radius which is related to the 
scalar curvature by the formula R = —2/a 2 . Furthermore, the volume form equals 
dV = a smh(g/a)dg A d<f>. The Hamiltonian for a free particle of mass m = 1/2 
takes the form 

rrO ( A 1 \ 1 9 _ u 1 

where Alb is the Laplace-Beltrami operator and g — det g^. We have set H = 1. 

The choice of a potential modeling the confinement is ambiguous. We nat- 
urally require that the potential takes the standard form of the quantum dot 
potential in the flat limit (a — > oo). This is to say that, in the limiting case, it 
becomes the potential of the isotropic harmonic oscillator Voo = j^ 2 g 2 . However, 
this condition clearly does not specify the potential uniquely. Having the freedom 
of choice let us discuss the following two possibilities: 

a) K(f) = ia 2 ^ 2 tanh 2 f, (2.1) 

b) U a (g) = \a 2 u, 2 sinh 2 £ (2.2) 

Potential V a is the same as that proposed in [7] for the classical harmonic 
oscillator on the Lobachevsky plane. With this choice, it has been demonstrated 
in [7] that the model is superintegrable, i.e., there exist three functionally inde- 
pendent constants of motion. Let us remark that this potential is bounded, and 
so it represents a bounded perturbation to the free Hamiltonian. On the other 
hand, the potential U a is unbounded. Moreover, as shown below, the stationary 
Schrodinger equation for this potential leads, after the partial wave decomposition, 
to the differential equation of spheroidal functions. The current paper concentrates 
exclusively on case b). 

The impurity is modeled by a (5-potential which is introduced with the aid 
of self-adjoint extensions and is determined by boundary conditions at the base 
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point. We restrict ourselves to the case when the impurity is located in the center 
of the dot (g = 0). Thus we start from the following symmetric operator: 

H = - GS + \ coth (!) I + h sinh " (!) a v sinh2 (f 



Dom(ff) = C£°((0,oo) x S 1 ) C L 2 f(0,oo) x S x ,a sinh(-) dpd0 



(2.3) 



2.2. Partial wave decomposition 

Substituting £ = cosh(g/a) we obtain 



d4> 2 



—(e-i)- 1 - 

4 ^ 7 4 



Bom(H) = C* °°((l,oo) x S 1 ) C L 2 ((l,oo) x S\a 2 d£ 



(2.4) 



Using the rotational symmetry which amounts to a Fourier transform in the vari- 
able <f>, H may be decomposed into a direct sum as follows 



H 



(J) H m , 



H„ 



d 



d 



4, ,2 



a u> 



i 



(e 2 



Dom(H m ) = Cg°(l,oo) C L 2 ((l,oo),d£). 
Note that i? m is a Sturm-Liouville operator. 

Proposition 2.1. _ff m is essentially s elf- adjoint for m ^ ; Ho deficiency indices 
(1,1). 

Proof. The operator is symmetric and semibounded, and so the deficiency 
indices are equal. If we set 

a 4 cj 2 1 
M = M, 461 = — , A = -z - - , 

then the eigenvalue equation 



_9 



(i 2 



1) 



d_ 



i 



(e 2 - 1) 



-0 = Zip 



(2.5) 



takes the standard form of the differential equation of spheroidal functions (|A.1[) . 
According to chapter 3.12, Satz 5 in 6J, for [i = m € No a fundamental system 
{yii J/ii} of solutions to equation (|2.5p exists such that 



yi(0 = (i-0 m/2 W-0, Vi(0) = i, 

wn(0 - (i - er m/2( p 2 (i - o + A m W (o bg (i - o, 
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where, for |£ — 1| < 2, ^i,*}^ are analytic functions in £, A, 8; and A m is a 
polynomial in A and 9 of total order to with respect to A and y/6; Aq = —1/2. 

Suppose that z € C\M. For to = 0, every solutions to (|2.5p is square integrable 
near 1; while for to =^ 0, y\ is the only one solution, up to a factor, which is square 
integrable in a neighborhood of 1. On the other hand, by a classical analysis due 
to Weyl, there exists exactly one linearly independent solution to (|2.5p which is 
square integrable in a neighborhood of oo, see Theorem XIII. 6. 14 in |S]. In the 
case of to = this obviously implies that the deficiency indices are (l,l).lfm^0 
then, by Theorem XIII. 2. 30 in [8], the operator H m is essentially self-adjoint. □ 



Define the maximal operator associated to the formal differential expression 

£ -!(«■- 4) +! T ! <«'- 1 '-i 

as follows 

T>om{H max ) = j/ G L 2 ((l,oc),d£) : /, /' G AC((l,oo)), 



- ^ - !) g J + — « - !)/ e L «>). d 
According to Theorem 8.22 in [9], H max = Ht. 

Proposition 2.2. Let k G (—00,00]. TTie operator Hq(k) defined by the formulae 
Dom(ff (K)) = {/ G Dom(i? ma:c ) : f 1 = nf } , H (n)f = H max f, 

where 

* : = - W A^rV lo^l-i)) - * : = A m + /0 log ( 2a ^ - ^ 

is a self-adjoint extension of H . There are no other self-adjoint extensions of H . 

Proof. The methods to treat 6 like potentials are now well established [TO]. Here 
we follow an approach described in [IT] , and we refer to this source also for the 
terminology and notations. Near the point £ = 1, each / G Dom(H max ) has the 
asymptotic behavior 

/(0 = /o^,l) + /i + o(l) as£^l+ 

where ,/o,/i G C and i* 1 (£,£') is the divergent part of the Green function for the 
Friedrichs extension of Hq. By formula (|2. 1 1|) which is derived below, F(£,l) = 
-l/(47ra 2 )log(2a 2 (£- 1)). Proposition 1.37 in [TTJ states that (C,ri,r 2 ), with 
IT/ = /o and r 2 / = /1, is a boundary triple for H max . 

According to theorem 1.12 in [TTJ, there is a one-to-one correspondence be- 
tween all self-adjoint linear relations k in C and all self-adjoint extensions of H a 
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given by k < — ► Hq(k) where Hq(k) is the restriction of H max to the domain of 
vectors / € DoTa(H max ) satisfying 

(Tif,T 2 f)eK. (2.6) 

Every self-adjoint relation in C is of the form k = Cv C C 2 for some v € R 2 , v ^ 0. 
If (with some abuse of notation) v = (1, k), k € M, then relation (|2.6[) means that 
fi = k/o. If w = (0, 1) then (|2.6p means that fo — which may be identified with 
the case k = oo. □ 



Remark. Let qo be the closure of the quadratic form associated to the semibounded 
symmetric operator Hq. Only the self-adjoint extension Hq(oo) has the property 
that all functions from its domain have no singularity at the point £ = 1 and 
belong to the form domain of qo- It follows that Hq(oo) is the Friedrichs extension 
of H (see, for example, Theorem X.23 in [T^] or Theorems 5.34 and 5.38 in [9]). 

2.3. The Green function 

Let us consider the Friedrichs extension of the operator H in L 2 ((1, oo) x S 1 , d£ dip) 
which was introduced in (|2.4p . The resulting self-adjoint operator is in fact the 
Hamiltonian for the impurity free case. The corresponding Green function Q z is 
the generalized kernel of the Hamiltonian, and it should obey the equation 



oo 

(H-z)g,(Z,<kt , ,tl/) = 8(Z-t , )6(<l> -<(/) = — J2 S(£-?)j m «-*"> 



m— — oo 



If we suppose Q z to be of the form 

^ oo 

e*(£,0;£',</>') = T- E G?(^'y m{ *-*'\ (2-7) 

m— — oo 

then, for all rn G Z, 

{H m -z)g?{t, a') = Ki-0- (2-8) 

Let us consider an arbitrary fixed and set 

a 4 uj 2 , 1 

n = to, AO = — , A = — z — -. 

4 4 

Then for all £ ^ equation (|2.8|) takes the standard form of the differential 
equation of spheroidal functions (jA.lj) . As one can see from (|A.8|) . the solution 
which is square integrable near infinity equals S'i™ 1 "' 3 '' (^, — a 4 w 2 /16). Furthermore, 
the solution which is square integrable near £ = 1 equals PsL™' (£, — a 4 CL> 2 /16) as 
one may verify with the aid of the asymptotic formula 

P ?(0 ~ 9m i^ + I)" 72 ase-l+, for tog N . 

2 m ' * ml 1 [y — m + 1) 
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We conclude that the mth partial Green function equals 

ei(f,f,_ (f-w.i m u w|W ) w r i6 ) Sv l e> ' i6 

(2.9) 

where the symbol W(Psu\Sv ) denotes the Wronskian, and £<, £> are re- 
spectively the smaller and the greater of £ and By the general Sturm-Liouvillc 
theory, the factor (£ 2 - l)W(Ps 1 " 11 , si" l|(3) ) is constant. Since Qf = Q- m decom- 
position (|2.7[) may be simplified, 

1 00 

= + - £^(£,£0 cos K4> -</»')]• (2-10) 

Z7T 7T * — ' 

m— 1 

2.4. The Krein Q-function 

The Krein Q-function plays a crucial role in the spectral analysis of impurities. It 
is defined at a point of the configuration space as the regularized Green function 
evaluated at this point. Here we deal with the impurity located in the center of 
the dot (£ = 1, (f> arbitrary), and so, by definition, 

Q{z) :=er»(l,0;l,0). 
Due to the rotational symmetry, 

&(0 := <t>\ 1, 0) = G z (t 0; 1, 0) = g,{t, 0; 1,0) = ^ G° z (t 1), 
and hence 

(#o-^(0=0, far£e(l,oo). 
Let us note that from the explicit formula (12. 9|) , one can deduce that the coefficients 
G™(£, 1) in the series in (|2.10|) vanish for m = 1,2,3,.... The solution to this 
equation is 

The constant of proportionality can be determined with the aid the following 
theorem which we reproduce from [13] . 

Theorem 2.3. Let d(x, y) denote the geodesic distance between points x,y of a 
two-dimensional manifold X of bounded geometry. Let 

n 

U S V(X) := (if : U+ := max([/, 0) e Lf o c (X), f/_ := max(-[/,0) e^L Pi (^)} 

i=l 

/or an arbitrary n G N and 2 < Pi < oo. Then the Green function Qjj of the 
Schrodinger operator Hjj = —Alb + U has the same on-diagonal singularity as 
that for the Laplace- Beltrami operator itself i.e., 

Gu{C,x,y) = -i- log 1 +Gu 9 ((]x,y) 
2?r d(x,y) 

where Q T f 9 is continuous on X x X . 
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Let us denote by and Q H (z) the Green function and the Krein Q-function 
for the Friedrichs extension of H, respectively. Since H = a 2 H and (H — z)Q z = 6, 
we have 

Qf<£, </>; e, </>') = a 2 g a 2 z (£, <j>; </>'), Q H (z) = a 2 Q(a 2 z). 
One may verify that 

logd(£>,0;0) = logg = log(aargcosh£) = ilog(2a 2 (£ - 1)) + 0(£ - 1) 
as g — > 0+ or, equivalently, £ — > 1+. Finally, for the divergent part 

F(t o == 0; e', - er 9 (e, 0; e', <£) = o ; e', 0) - er 9 (^ o ; 0) 

of the Green function Q z we obtain the expression 

^,1) = -^ log(2a 2 (e-l)). (2.11) 

From the above discussion, it follows that the Krein Q-function depends on the 
coefficients a, (3 in the asymptotic expansion 

5 " (3) (£<-^r) =alog(£-l)+/3 + o(l) as^l+, (2.12) 

and equals 

v ; 47ra 2 a 47ra 2 

To determine a, /? we need relation (|A.10p for the radial spheroidal function 
of the third kind. For v and v + 1/2 being non-integer, formula (|A.12[1 implies that 

S^H^O) = ^M e -.vK^o (Ws o^ iK|fl)) 

x (2-14) 

7T 

Applying the symmetry relation (|A.5[) for expansion coefficients, we derive that 

oo 

Q*°,,-i(e,0) = X; (-l) r a°_ u _ 1>r (9)Q°_ v _ 1+2r (0 

r— — oo 
oo 



r— — oo 



Using the asymptotic formulae (see [5]) 

Ql(0 = -\ log ^ + *(i) - + 1) + o((£ - i) io g (e - 1)) 
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the series expansion in (|A.11[) and formulae (|2 . 14[> . we deduce that, as £ — * 1+, 



\{0)- x U lo s 



*(1) + Trcot(^Tr) + ^s v {0) 



7T 

4m- 1 Qbg^ >imi*: 



where the coefficients are introduced in (|A.7p , 

oo 

*s v (0):= (-l)'<,W*(^ + l + 2r) ! 

r— — oo 

and where we have made use of the following relation for the digamma function: 

&(—z) = + 1) + 7TCOt(7Tz). 

We conclude that 



S° (3) (£, B) ~ a log(C - 1) + p + O ((£ - 1) log(C - 1)) as £ - 1+, 



where 



a 



i tan(^7r) 



(e™ " e^ 2 ^ 3 / 2 ^^)) , 



2tT S °(0) 

/? = a (- log 2 - 2* (1) + 2¥a„(0)s°(0)) + c 
The substitution for a, /3 into (|2.13p yields 
1 



Q(*)=- 



47ra 2 



- log 2 - 2*(1) + 2 *s J - 



16 



1G 



2a 2 tan(^7r) I 
where i/ is chosen so that 



^7r(3i>+3/2) 16 I 



A° - 



- 1 



log(2a 2 
Ana 2 



16 



(2.15) 



(2.16) 



For v — n being an integer, we can immediately use the known asymptotic 
formulae for spheroidal functions (see Section 16.12 in [5]) which yield 



" +ir law } ' 



2r>— n 

as £ -> 1+. Here, h = 1, /ife = 1/1 + 1/2 + . . . + 1/k. By (|2~T5|) . one can calculate 
the Q-frmctkm in this case, too. 
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2.5. The spectrum of a quantum dot with impurity 

The Green function of the Hamiltonian describing a quantum dot with impurity 
is given by the Krein resolvent formula 

0f« (£, <j>; 0') = G? & 0; - ggj) (f, 0; 1, O)0f (1, 0; 0) 

(recall that, due to the rotational symmetry, G% (£, <^>; 1, 0) = G^(£,, 0; 1, 0)). The 
parameter x := a 2 '* € (—00,00] determines the corresponding self-adjoint exten- 
sion H(x) of H. In the physical interpretation, this parameter is related to the 
strength of the 5 interaction. Recall that the value \ — 00 corresponds to the 
Friedrichs extension of H representing the case with no impurity. This fact is also 
apparent from the Krein resolvent formula. 

The unperturbed Hamiltonian iJ(oo) describes a harmonic oscillator on the 
Lobachevsky plane. As is well known (see, for example, [14] ). for the confinement 
potential tends to infinity as g — > 00, the resolvent of H (00) is compact, and the 
spectrum of H(oo) is discrete and semibounded. The eigenvalues of H(oo) are 
solutions of a scalar equation whose introduction also relies heavily on the theory 
of spheroidal functions. We are sceptic about the possibility of deriving an explicit 
formula for the eigenvalues. But the equation turned out to be convenient enough 
to allow for numerical solutions. A more detailed discussion jointly with a basic 
numerical analysis is provided in a separate paper |15j . 

A similar observation about the basic spectral properties (discreteness and 
semiboundedness) is also true for the operators H(\) for any x£K since, by the 
Krein resolvent formula, the resolvents for H(x) and H(oo) differ by a rank one 
operator. Moreover, the multiplicities of eigenvalues of H(x) an d H(oo) may differ 
at most by ±1 (see j9j Section 8.3]). 

A more detailed and rather general analysis which is given in [1] can be carried 
over to our case almost literally. Denote by a the set of poles of the function Q H (z) 
depending on the spectral parameter z. Note that a is a subset of spec(i? (00)). 
Consider the equation 

Q H (z) = x- (2.17) 

Theorem 2.4. The spectrum of H(x) * s discrete and consists of four nonintersecting 
parts S±, S2, S3, S4 described as follows: 

1. Si is the set of all solutions to equation \2.17\ ) which do not belong to the 
spectrum of H (00). The multiplicity of all these eigenvalues in the spectrum 
of H(x) equals 1. 

2. 52 is the set of all A € a that are multiple eigenvalues of H (00). If the mul- 
tiplicity of such an eigenvalue A in spec(i/(oo)) equals k then its multiplicity 
in the spectrum of H(x) equals k — 1. 

3. S3 consists of all A G spcc(i?(oo))\cr that are not solutions to equation \2.11^ . 
the multiplicities of such an eigenvalue A in spec(i/(oo)) and spec(i?(%)) are 
equal. 
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4. 64 consists of all A € spec(7?(oo)) \ a that are solutions to equation £2.17)) . 
If the multiplicity of such an eigenvalue A in spec(_ff (00)) equals k then its 
multiplicity in the spectrum of H{\) equals k + 1. 

Hence the eigenvalues of H(x), x G K, different from those of the unperturbed 
Hamiltonian H(oo) are solutions to (|2.17p . As far as we see it, this equation can 
be solved only numerically. We have postponed a systematic numerical analysis of 
equation (|2.17p to a subsequent work. Note that the Krein Q-function (|2.15|) is in 
fact a function of v, and hence dependence (|2.16p of the spectral parameter z on 
v is fundamental. In this context, it is quite useful to know for which values of v 
the spectral parameter z is real. A partial answer is given by Proposition IA. l[ 



3. Conclusion 

We have proposed a Hamiltonian describing a quantum dot in the Lobachevsky 
plane to which we added an impurity modeled by a S potential. Formulas for the 
corresponding Q- and Green functions have been derived. Further analysis of the 
energy spectrum may be accomplished for some concrete values of the involved 
parameters (by which we mean the curvature a and the oscillator frequency u>) 
with the aid of numerical methods. 



Appendix: Spheroidal functions 

Here we follow the source [5]. Spheroidal functions are solutions to the equation 

(i - a§ - 2^ + [a + 40a - e) - - e 2 )- 1 ] * = 0, (a.i) 

where all parameters are in general complex numbers. There are two solutions 
that behave like times a single- valued function and times a single- valued 

function at 00. The exponent v is a function of A, 9, fi, and is called the charac- 
teristic exponent. Usually, it is more convenient to regard A as a function of ^, fi 
and 9. We shall write A = A£(0). If v or fi is an integer we denote it by n or m, 
respectively. The functions A£J(#) obey the symmetry relations 

K{9) = = AVr(0) = AliUP). (A.2) 

A first group of solutions (radial spheroidal functions) is obtained as expan- 
sions in series of Bessel functions, 

00 

sFHw) = C 1 - r 2 )- M/a *s(0) E <A%rW 1/2 o, (a.3) 
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j = 1,2,3, 4, where the factors s^(9) are determined below and 



2^ --«+i/2V.s;> v-v vs> -y 2C -",,+1/2 

The coefficients a^ r (#) (denoted only a r for the sake of simplicity) satisfy a three 
term recurrence relation 

(v + 2r - fi)(v + 2r - fj,- 1) „ [y + 2r + fi + 2){v + 2r + fi + 1) 

■ Ud r —\ H — ~ /rtX / — /rtX C/G r -|-1 



(z/ + 2r-3/2)(^ + 2r- 1/2) 



A£(0)-(^ + 2r)(^ + 2r + l) + 



(i/ + 2r + 3/2)(i/ + 2r + 5/2) 



(u + 2r)(v + 2r + 1) + -1 
(i/ + 2r-l/2)(^ + 2r + 3/2) 



20 



a,. = 0. 
(A.4) 



Here and in what follows we assume that v+1/2 is not an integer (to our knowledge, 
the omitted case is not yet fully investigated). 

The coefficients r (9) may be chosen such that 



<o(9) = ^-i,oW = ^gW, 



and so (see (IA.2I0 



[1/ + /i + lj2r 



(A.5) 



where (a) r := a(a + l)(a + 2) . . . (a + r — 1) = T(a + r)/T(a), (a)o := 1. Equation 
(|A.4[) leads to a convergent infinite continued fraction and this way one can prove 
that 

CL r T dr 9 

lim = lim = — . (A. 6) 

r^oo a r _i r— »— 00 JZf-i-i 4 

From (|A.6|) and the asymptotic formulae for Bessel functions, it follows that (I A.3|) 
converges if |£| > 1. 
If we set in (|A.3|) 

-1 



(A.7) 



then 



SgU)({,0) ~ V^(26< 1/2 0, for I &rg{9 1/2 0\ as £ ^ 00. 

We have the asymptotic forms, valid as £ — * 00, 

<?M 3 )(£,0) = i^-V2 r 1^(2^-^/2^/2) [1 + (| f |-l)] j 

for - 7T < arg(6» 1/2 < 2tt, 
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and 

SfM(£,0) = 1 tf -l/2 € _l e -i(2flV>e_ W2 _ w/2)[1 + (|^|-l)] j 

2 (A. 9) 

for - 2tt < arg((9 1/2 £) < tt. 

The radial spheroidal functions satisfy the relation 

S£( 3 ) = 1 —— (s^L + i e~ inv S^) . (A.10) 

i cos(vn) V / 

The radial spheroidal functions are especially useful for large £; the larger is 
the £ the better is the convergence of the expansion. To obtain solutions useful near 
±1, and even on the segment (—1, 1), one uses expansions in series in Legendre 
functions, 

r= ~°° (A.H) 

r— — cjo 

These solutions are called the angular spheroidal functions and are related to the 
radial spheroidal functions by the following formulae: 

Sp\t, 9) = tt- 1 sin[(^ - ^e-^^K^Qs^^, 9), 

where K£{6) can be expressed as a series in coefficients a$ r (9), and sometimes it 
is called the joining factor. In more detail, for any k € Z it holds true that 

K^(9) = \ (JJ ^ r(l + v — [i + 2k) e^+^stie) 

(~ 1 ) r <r(«) 



(fc-r)!r(^ + fc + r + 3/2) 



r— 



(r-Jfe)!r(l/2-i/-Jfe-r) 



Proposition A.l. Let i/,6 G M and set /u = 0. TTien A® (9) £ 
Proof. To simplify the notation we denote, in (|A.4[) . 

{v + 2r + [i + 2)(y + 2r + fi + 1) 



US T &l TO/'MfT ^' 



(^ + 2r + 3/2)(^ + 2r + 5/2) 

(^ + 2r)(^ + 2r + l) + ^ 

(y + 2r - fi){v + 2r - fi - 1) 
+ 2r-3/2)(z/ + 2r- 1/2) 
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The resulting formula may be written in the matrix form, 



/ • 






\ 


{ '■■ ) 




f : \ 




7-i P-i ot-i 






ct-i 




ct-i 




7o A) 


Of) 




a 


= -A 






7i 






ai 




ai 


V 






■■) 


V : ) 




V : / 



(A.13) 



where we have omitted the fixed indices. 
As one can see, 



„,„ (f) , _ v + 2r + 5/2 
^^-^ 2^1/2^ W 

and so 

^ + 2r - 3/2 a °- l(g)a °^-i (g) +^" / ( g Xr( g ) + «?'' / ( g )<r + i(g) = - A°(g)o° r (g). 
Substitution a° r = L r (v)aP vr , where L r (y) are non-zero constants, yields 
^ + 2r+l/2 o ^L r _i(i>) i Q0>l//m . , m , ^.o.r/^-o 



f^i^Xr-xW^pT^ + #'"(0)5° r (0) + a^(0)a°, r+1 (0) : 

Li r \y) 



L r (v) 



v + 2r-'3/2 
= -A°(0)a° r (0). 

We require the matrix in (|A.13|) to be symmetric in the new coordinates {a r }. 
This implies that 

v + 2r + 1/2 L r - X (v) _ L r (v) 
v + 2r-3/2 L r (v) ~ L r -x(v)' 



For r $ {-v/2 - 1/4, -v/2 + 3/4), the solution is L r (v) = sj\v + 2r + 1/2]. For 
ro = re (—v/2— 1/4, — v/2 + 3/4), there is no real solution and so we set L rQ (v) — 
yJzT-j- 2ro + 1/2| and make another transformation of coordinates: 

- _f-a r for r = r - (2fc - 1), k £ N, 
1 a r for all other r. 

Relation (|A.13j) can be viewed as an eigenvalue equation with a symmetric 
matrix in the coordinate system {cifc}, hence A°,(#) must be real. □ 
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